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1 . INTRODUCTION 

The  objective  of  this  program  was  to  develop  and  demon- 
strate an  energy  source  having  up  to  several  hundred  times 
the  energy  density  of  high  explosives.  Such  a high  energy 
density  source  would  provide  an  Important  addition  to  nuclear 
burst  simulation  technology.  For  example,  the  source  energy 
density  Is  of  primary  Importance  In  simulating  the  high 
temperature  ablating  gas  flow  for  a nuclear  source  In  a 
trench.  Situations  In  which  the  size  and  mass  of  the  source 
are  Important  Include  the  alrblast  coupling  from  a shallow 
penetrating  nuclear  burst  and  the  bubble  dynamics  of  an 
underwater  nuclear  burst.  Because  of  the  large  difference  In 
energy  density  between  a nuclear  device  and  chemical  explosive 
(approximately  k orders  of  magnitude)  It  Is  Inadequate  to 
simulate  these  events  with  a chemical  explosive. 

The  program  was  specifically  concerned  with  the  experi- 
mental development  of  an  explosive  magnetohydrodynamlc  (MHD) 
generator  to  concentrate  and  convert  the  energy  of  chemical 
explosives  Into  electrical  form.  The  electrical  energy  can 
be  transmitted  to  an  exploding  wire  achieving  an  energy 
density  at  least  100  times  greater  than  the  chemical 
explosive.  While  It  Is  recognized  that  the  anticipated 
energy  density  Is  still  not  a direct  simulation  of  a nuclear 
source,  the  expected  difference  In  energy  density  would 
provide  an  opportunity  to  perform  alrblast  coupling  experl- 
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ments  at  shallower  burial  depths  than  now  possible  as  well 
as  provide  information  on  the  dependence  of  ground  motion 
and  alrblast  phenomena  on  source  energy  density.  The  unique 
feature  of  using  an  explosively  driven  MHD  generator  to 
provide  the  electrical  energy  is  that  the  device  can  be 
fabricated  in  any  size  scale  requiring  up  to  several  hundred 
tons  of  explosive  to  simulate  the  energy  equivalent  of  up  to 
several  tons  of  chemical  explosive.  Alternate  high  power 
energy  sources  such  as  capacitor  banks  are  limited  in  total 
stored  energy  to  several  megajoules  and  are  prohibitively 
expensive  for  a large  amount  of  stored  energy.  For  example 
the  AURORA  capacitor  bank  stores  about  5 megajoules,  which 
at  a 2 megajoule  per  pound  TNT  yield  equivalence  contains 
the  energy  of  only  2 1/2  pounds  of  explosive.  A goal  of 
this  program  was  to  deliver  20  megajoules  of  energy  to  an 
exploding  wire. 

A schematic  of  the  exploding  wire  simulator  test  con- 
cept is  shown  In  Figure  1.1.  It  consists  of  an  explosive 
generator,  a transmission  line,  and  an  exploding  wire  simu- 
lator. The  explosive  generator  Is  located  a considerable 
distance  from  the  simulator  source  so  that  blast  effects 
from  the  generator  do  not  affect  simulator  measurements. 

The  required  size  of  the  generator,  which  directly  affects 
test  cost.  Is  governed  by  overall  energy  coupling  efficiency 
between  chemical  energy  In  the  explosive  and  delivered  energy 


FIGURE  l.l  SCHEMATIC  OF  EXPLODING  WIRE  SIMULATOR  TEST  CONCEPT 


in  the  exploding  wire  simulator.  This  coupling  may  be 
expressed  as  a ratio  of  yield  (pounds  of  TNT  equivalent)  to 
pounds  of  explosive  required. 


The  actual  electrical  energy  deposited  in  the  exploding 
wire  can  be  determined  using  measurements  of  current  and 
voltage.  This  energy  represents  the  total  energy  available 
for  air  blast,  ground  motion,  vaporization  of  the  surrounding 
media  and  other  source  phenomena.  The  equivalent  yield  of 
the  exploding  wire  source  Is  typically  of  interest  In  nuclear 
effects  studies  and  Is  best  determined  by  fitting  obtained 
data  (alrblast,  ground  motion,  etc.)  with  scaled  nuclear 
data.  The  equivalent  yield  of  the  source  Is  the  scaled 
nuclear  yield  which  best  fits  the  data. 

Overall  coupling  between  the  explosive  In  the  generator 
and  the  equivalent  yield  of  the  exploding  wire  source  can 
only  be  determined  experimentally,  as  it  depends  on  many 
Interrelated  factors  such  as  explosive  driver  efficiency, 

MHD  generator  efficiency,  coupling  between  the  generator 
and  the  load,  losses  In  the  transmission  line,  and  so  on. 
Rough  estimates,  supported  by  the  discussions  below.  Indicate 
coupling  efficiencies  In  the  range  5 pounds  per  pound  to  50 

I 

pounds  per  pound.  The  objective  of  the  present  program  was 
restricted  to  demonstrating  the  feasibility  of  the  explosive 

! MHD  generator.  Experimental  Investigation  of  transmission 

I 

line  and  exploding  wire  coupling  will  be  the  objective  of 
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' later  programs.  A specific  program  goal  was  to  deliver  20 

megajoules  of  energy  to  an  exploding  wire. 


2.  TECHNICAL  DISCUSSION 
2.1  Explosive  MHD  Generator 

The  explosive  VHD  generator  used  In  this  program  com- 
prises two  elements:  an  explosive  driver,  which  transfers 
and  concentrates  chemical  energy  of  a high  explosive  Into 
a highly  Ionized  and  electrically  conductive  Intermediary 
gas;  and  an  MHD  channel,  which  converts  the  energy  of  the 
conductive  gas  Into  electrical  energy. 

A schem.atlc  of  explosive  driver  operation  Is  shown  in 
Figure  2.1.  It  consists  of  a tube  filled  with  explosive 
enclosed  in  a pressure  tube  contjining  driver  gas.  When 
the  explosive  Is  detonated  at  one  end,  the  expansion  tube 
expands  and  Impacts  the  pressure  tube;  this  impact  produces 
a moving  gas  seal  which  drives  the  gas  forward  to  form  a 
high  pressure  shock  wave.  Energy  Is  transferred  from,  the 
explosive  to  the  driver  gas  because  the  explosive  gases  must 
do  work  against  the  driver  shock,  pressure  in  order  to  expand. 
The  explosive  driver  is  potentially  very  efficient;  maximum 
efficiency  would  be  obtained  when  the  detonation  gases  do  so 
much  work  that  the  expansion  tube  Impact  velocity  goes  to 
zero.  Practical  limits  are  Imposed  by  the  requirement  of 
maintaining  a gas  seal,  which  depends  In  a complex  fashion 
on  Impact  velocity,  tube  materials,  and  dynamic  boundary 
layer  phenomena  In  the  gas.  Reasonable  extrapolations  Indi- 
cate that  50%  efficiency  can  be  expected;  achieved  efficiency 
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to  date  (no  experimental  program  dedicated  to  optimizing 
efficiency  has  ever  been  carried  out)  Is  2^%.  The  conser- 
vative driver  selected  for  the  current  program  delivers 
about  10%  efficiency. 

Gasdynamlc  conditions  achieved  by  an  explosive  driver 

are  uniquely  suited  to  MHD  electric  power  generation.  The 

energy  density  of  the  gas  (typically  about  72  kilojoules 

per  gram  total  energy)  Is  sufficiently  high  that  the  gas  Is 

completely  Ionized  without  the  requirement  for  seed  material 

or  other  exotic  techniques  In  common  use  with  combustion 

driven  MHD.  Conductivity  ranges  from  about  300  mho/cm  for 

argon  to  about  800  mho/cm  for  xenon,  or  about  the  conductivity 

of  graphite.  Gas  velocities  are  typically  on  the  order  of 

1 cm/  /xsec,  and  mechanical  power  flow  on  the  order  of  10^^ 

2 

watts/cm  . The  mechanical  flow  conditions.  Including  shock 
velocity,  power  flow  rate,  etc.,  have  been  verified  many 
times  In  numerous  programs  using  explosive  drivers.  The 
electrical  properties  are  being  Investigated  In  the  present 
program. 

A simplified  schematic  of  the  MHD  generator  concept 
used  In  this  program  Is  Illustrated  In  Figure  2.2.  The 
device  operates  as  a Faraday  crossed- field  generator,  with  a 
highly  conductive  hot  gas  flowing  normal  to  a magnetic  field. 
According  to  Maxwell's  equations,  an  electrical  field  Is 
Induced  In  the  gas  In  a direction  mutually  orthogonal  to 
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FIGURE  2.2  SELF-EXCITED  MHD  GENERATOR 


the  gas  velocity  and  the  magnetic  field.  Electrodes  exposed 
to  the  gas  flow  and  aligned  with  the  direction  of  the  electric 
field  can  be  used  to  measure  the  voltage  Induced  across  the 
gas  channel  and  can  be  used  to  energize  an  external  circuit. 

The  electrical  energy  delivered  to  an  external  load 

depends  upon  the  complex  Interaction  between  the  magnetic 

field  and  the  conductive  gas,  as  well  as  the  electrical 

v<*?iaracterlstlcs  of  the  load.  For  typical  gas  conditions 

produced  in  an  explosive  driver,  the  total  thermal  power 

per  unit  cross-sectional  area  of  the  flowing  gas  Is  1.1  x 
10  2 

10  watts/cm  . In  order  to  convert  a reasonable  portion 
of  gas  energy  Into  electrical  energy.  It  Is  necessary  to 
provide  a very  strong  magnetic  field  In  the  generator,  on 
the  order  of  10  Tesla.  This  field  strength  Is  well  beyond 
the  capabilities  of  permanent  magnets  ( ~ 0.1  Tesla)  and 
simple  electromagnets  (~1  Tesla),  but  can  be  achieved  In 
a generator  by  the  self-excltatlon  process  depicted  In 
Figure  2.2.  A portion  of  the  electrical  output  of  the  gen- 
erator Is  fed  through  a coll  located  along  the  gas  channel 
so  as  to  Increase  the  magnetic  field  strength  in  the  region 
of  the  electrodes.  The  Increased  magnetic  field  causes  the 
generator  voltage  to  Increase,  thereby  Increasing  the 
electrical  power  produced  by  the  generator.  The  process  of 
utilizing  a portion  of  the  generator  output  to  Increase  the 
field  strength  and  thereby  Increase  generator  output  Is 
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called  self-excitation  and  is  used  in  the  operation  of  all 
conunercial  power  generators. 

The  primary  technical  problem  addressed  in  this  program 
was  accomplishing  the  self-excitation  process  on  a micro- 
second time  scale.  The  successful  operation  of  a high-power 
self-excited  MHD  generator  for  driving  an  exploding  wire 
requires  a careful  design  and  integration  of  the  explosive 
and  gasdynamlc  system  with  the  electrodes,  self-excitation 
colls  and  load  characteristics.  A combined  analytical  and 
experimental  effort  was  conducted  to  match  the  electrical 
parameters  of  the  generator  with-  gasdynamlc  flow  parameters 
in  such  manner  as  to  maximize  the  electrical  energy  delivered 
to  an  exploding  wire.  Significant  effort  was  expended  to 
determine  the  electrical  and  flow  properties  of  the  gas  and 
to  understand  their  role  in  the  rapid  self-excitation  MHD 


interaction. 
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3.  EXPERIMENTAL  STUDY  OF  EXPLOSIVE  MKD  GENERATOR 

3.1  Gasdynamlc  Design  and  Performance 

3.1.1  Explosive  Driver 

The  purpose  of  the  explosive  driver  Is  to  transfer 
explosive  energy  to  a selected  driver  gas,  shock-compressing 
It  to  a hlgh-temperature  electrically  conducting  plasma 
state.  Efficiency  of  the  energy  transfer  process  Is  an 
important  factor  for  large  tests,  as  it  affects  the  amount 
of  explosive  required  per  unit  of  simulator  source  yield, 
but  Is  less  Important  for  smaller  tests.  Our  approach  In 
this  program  was  to  select  a relatively  simple  non-optimum 
driver  configuration  with  limited  efficiency  and  concentrate 
effort  on  the  MHD  process. 

Prior  experience  has  indicated  that  a driver  comprising 
an  aluminum  expansion  tube  and  a steel  outer  tube  Is  moder- 
ately effective.  Aluminum  can  expand  substantially  without 
' rupture,  and  aluminum-steel  Impact  pressures  are  sufficiently 

i 

j high  to  overcome  boundary  layer  stagnation  pressures  in  the 

r driver  gas,  which  may  be  many  times  higher  than  design  shock 

I pressure.  For  a given  design,  maximum  efficiency  Is  obtained 

r 

by  operating  the  driver  at  as  high  a pressure  as  possible; 
the  limit  Is  governed  by  loss  of  gas  seal  at  high  shock 
pressures . 

Driver  configuration  for  the  program  Is  shown  In  Figure 
3.1,  and  design  parameters  based  on  the  assumption  of  minimal 


FIGURE  3.1  SCHEMATIC  OF  EXPLOSIVE  DRIVER 


gas  leakage  are  shown  In  Table  3.1.  Data  from  Shot  120-1, 
which  had  a 10-kbar  dealg'n  pressure.  Indicate  excessive  gas 
leakage  on  the  order  of  0-501.  Rather  than  Investigate  the 
specific  reason  for  this  leakage  the  design  shock  pressure 
was  reduced  from  10  kbar  to  6 kbar  in  subsequent  shots. 

This  modification  reduced  the  leakage,  although  it  also 
reduced  design  efficiency.  Gas  leakage  was  determined  by 
comparing  the  actual  shock  location  with  ideal  shock  location 
based  on  Hugonlot  shock  relationships  and  known  geometric 
configurations . 

With  minimal  gas  leakage,  the  drivers  used  in  Shots 
120-2  to  120-4  delivered  about  10?  efficiency,  or  10  pounds 
of  explosive  yield  required  per  pound  of  equivalent  plasma 
energy  yield.  The  Hugonlot  shock  Jump  conditions  determine 
the  plasma  state,  and  since  measured  shock  velocity  in  the 
driver  was  very  close  to  ideal,  the  plasma  conditions  were 
close  to  calculated  values. 

3.1.2  Transition  Section  and  Channel 

Driver  plasma  is  produced  in  an  annular  shape,  and  is 
delivered  to  the  cylindrical  MHD  channel  by  means  of  a 
tapered  transition  section.  The  taper  also  serves  to 
Increase  pulse  duration  by  decreasing  cross-sectional  area. 
Shots  120-1  to  120-3  were  designed  to  have  an  area  reduction 
from  driver  annulus  to  channel  of  3.3  to  1,  and  Shot  120-4, 


In  an  attempt  to  increase  pulse  duration  even  more  to  com- 
pensate for  low  plasma  conductivity,  had  an  area  reduction 
of  7.2  to  1. 
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I The  driver  plasma  pulse  was  further  processed  by  expan- 

sion into  the  channel  prior  to  interaction  with  the  MHD 
generator.  Expansion  serves  the  dual  purpose  of  increasing 
pulse  duration  and  increasing  velocity  of  the  leading  edge 
of  the  pulse,  both  of  which  are  important  for  self-excited 
operation.  Shots  120-1  to  120-3  were  fired  with  1-atm 
flushed  argon  Initially  in  the  channel,  and  120-4  was  fired 

with  0.25-atm  air  to  Increase  peak  velocity  still  further. 

S 

Initial  performance  predictions  of  the  transition 
section  and  channel  were  preformed  by  using  Isentroplc 
nozzle  relationships  for  the  channel.  More  detailed  analysis 
was  made  using  a one-dimensional  finite-difference  code  with 
variable  area  option.  Typical  computed  flow  speed  histories 
at  several  channel  locations  are  presented  in  Figure  3-2. 

3.1.3  Performance  Data 

Performance  of  driver  and  transition  section  gasdynamlcs 
were  established  experimentally  by  analyzing  tlme-of-arrlval 
data  on  position-time  plots  of  the  shot,  as  typified  by 
Figures  3.3  and  3.4.  Time  of  arrival  was  established  by 
self-shorting  cap  pins  in  the  driver  and  by  various  channel 
diagnostics.  As  noted  above,  the  driver  shock  velocity 
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Flow  Speed  (cm//Asec) 


Time  (/xsec) 


FIGURE  3.2  CALCULATED  FLOW  SPEED  VERSUS  TIME  IN  THE  MHD  CHANNEL  AT 
SEVERAL  AXIAL  DISTANCES  FROM  THE  DIAPHRAGM  FOR  SHOT  120-4 
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Channel  Shock 
1.34  cm//isec 


Detonation 
0.85  cm/^sec 


FIGURE  3.4  POSITION-TIME  DIAGRAM  FOR  SHOT  120-4 


was  close  to  the  calculated  values.  Shock  propagation  through 
the  transition  section  Into  the  channel  was  found  to  be 


slower  than  the  finite  difference  code  calculations,  probably 
because  the  one-dimensional  modeling  Inadequately  represented 
the  two-dimensional  Interactions  in  the  transition  section. 
Although  this  source  of  inaccuracy  should  eventually  be 
eliminated,  it  caused  only  minor  inconvenience  in  the  present 
program,  as  measured  velocities  were  used  to  estimate  MHD 
performance . 

Velocity  history  in  the  channel  was  obtained  with  a 
magnetic  velocity  gage  (an  open-circuit  MHD  generator).  A 
typical  velocity  gage  record  is  shown  in  Figure  3.5.  The 
plasma  pulse  is  smooth,  with  reasonable  duration,  and  appears 
suitable  for  MHD  power  generation.  Shot  120-^  was  an 
Increase  of  about  50X  in  size  scale  to  lengthen  pulse 
duration.  The  data  verify  that  the  driver  gasdynaralcs 
performed  equally  well  at  the  larger  scale. 

3.2  Generator  Design  and  Performance 

The  principle  technical  problem  foreseen  in  this  pro- 
gram was  rapid  self-excitation  of  the  MHD  generator.  In 
order  to  convert  a reasonable  portion  of  gas  energy  into 
electrical  energy,  it  is  necessary  to  provide  a very  strong 
magnetic  field,  on  the  order  of  10-50  tesla.  Although  it 
is  possible  to  provide  this  field  with  external  energization 
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(e.g.  a capacitor  bank  or  another  explosive  MHD  generator), 
the  practical  requirements  of  large  scale  field  tests  strongly 
favor  self-excltatlon.  With  self-excltatlon  an  Initial 
field  provided  by  permanent  magnets  or  an  external  capacitor 
bank  can  be  amplified  to  the  desired  operating  field.  All 
commercial  rotary  generators  operate  In  this  manner. 

3.2.1  Single-State  Analysis 

The  technical  problem  areas  In  devising  an  effective 
self-excltatlon  circuit  are  most  reaclly  Introduced  by 
considering  a simple  one-stage  self-excited  generator,  as 
shown  In  Figure  2.2  The  generator  voltage  Is  given  by 

V=Bub-Rpl  (3.1) 

where  B Is  the  magnetic  field,  u Is  the  plasma  velocity,  j 

I 

b Is  the  Interelectrode  distance,  Rp  Is  the  plasma  resls-  j 

i 

tance,  and  1 Is  the  generator  current,  all  In  consistent  f 

k 

MKS  units.  An  Initial  magnetic  field  B^  Is  provided,  and  | 

any  additional  field  Is  proportional  to  the  field  current.  1' 

Equation  (3.1)  can  therefore  be  written 

V = Vo  + (Z-Rp)I  (3.2) 

where  = B^ub  . The  parameter  Z Is  the  Incremental  ratio  j 

i 

: j 
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2U 


of  open  circuit  voltage  to  field  coll  current,  and  represents 
a feedback  resistance  from  output  to  Input. 


Applying  the  generator  to  an  Inductive  load,  we  have 


V=L 


dt 


+ RlI 


(3.3) 


or 


L-^  = (Z-R)l  + Vo  (3.4) 

where  R = Rp  + Rq  total  circuit  resistance  and  l Is 

the  total  circuit  Inductance. 

The  generator  self-excltes  If  Z>R  , l.e.  the  feedback 
resistance  exceeds  the  circuit  resistance.  The  solution  of 
Equation  (3.4)  Is 


1 = 


~(Z-R)t 

L 


(3.5) 


= — = exp 


~(Z-R) 

L 


(3.6) 


The  self-excltatlon  process  will  continue  until  the  power 
drawn  Is  a significant  fraction  of  the  thermal  power  of  the 
plasma,  at  which  point  the  velocity  (and  hence  Z)  will 
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decrease.  This  slowdown,  usually  associated  with  shock 
waves  In  the  flow.  Is  known  as  choking.  It  has  been  shown  In 
large  non-equlllbrlum  MHD  generators  that  a power  extraction 
rate  of  approximately  20%  can  be  achieved  without  serious 
choking  effects  (Reference  3*1)  • This  same  extraction  rate 
was  used  as  an  upper  limit  In  our  performance  calculations, 
recognizing  that  the  very  different  plasma  conditions  In  a 
high  density  equilibrium  MHD  generator  may  alter  this  assumed 
rate  limit. 

The  principal  technical  problem  Is  to  obtain  a combi- 
nation of  Vq  , Z , R , and  l to  achieve  choke  conditions 
before  the  plasma  pulse  passes  through  the  generator, 
depends  on  the  Initial  field,  Z depends  on  the  product  of 
field  coll  B/1  ratio  and  plasma  velocity,  R depends  princi- 
pally on  plasma  conductivity,  and  L Is  a field  coll  parameter. 
Circuit  design  concentrates  on  field  coll  B/I  and  L , and 
gasdynamlc  design  concentrates  on  plasma  pulse  duration, 
velocity,  and  conductivity. 

Considerable  engineering  effort  was  expended  designing 
the  field  coll  assembly  to  minimize  L and  maximize  B/I  . 

Pinal  measured  parameters  for  the  2-lnch  channel  size  are 

-K 

L ■ 70  nanohenry,  B/l  = 1.2  x 10  tesla/amp.  At  a velocity 
of  1 cm/mlcrosecond  this  value  of  B/l  corresponds  to  Z ■ 5 
mllllohm.  Coll  resistance  Is  a small  fraction  of  a mllllohm, 
and  Is  negligible  compared  to  Z . If  plasma  resistance  were 
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also  negligible,  the  field  doubling  time  would  be 


T = — In2  = 9.7/iSec 


(3.7) 


Plasma  resistance,  if  it  is  too  large,  has  a major 
effect  on  field  amplification  for  a given  plasma  pulse 
duration.  Resistance  is  related  to  plasma  conductivity  by 
a measured  geometrical  factor:  R = .28/<r  , where  <r  is 
measured  in  mho/cm.  The  theoretical  value  of  <r  is  about 
^<00  mho/cm,  corresponding  to  R * 0.7  mllllohms,  low  enough 
not  to  be  a factor.  At  o-  * 50  mho/cm,  R * 5.6  mllllohms 
and  the  plasma  resistance  would  overcome  the  feedback 
resistance  and  self-excitation  would  not  occur.  Inter- 
mediate values  lead  to  self-excitation  of  the  fields  but  at 
a longer  time  scale.  At  100  mho/cm,  for  example,  R ■ 2.8 
mllllohm  and 


T = In2  = 32/Asec  (3.8) 

£“K 

For  a given  pulse  duration  field  amplification  depends 
exponentially  on  the  doubling  time,  which  is  a function  of 
plasma  conductivity.  Calculated  field  amplification  for  a 
100  microsecond  pulse  duration  is  shown  in  Figure  3.6.  With 
conductivities  anywhere  near  theoretical,  there  is  more  than 
enough  field  amplification  to  ensure  choke  conditions  even 
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from  a small  Initial  field.  Low  conductivities  on  the  order 
of  100  mho/cm  lead  to  greatly  reduced  amplification,  and  the 
actual  value  Is  highly  sensitive  to  achieved  conductivity. 

Although  analysis  of  multi-stage  circuits  with  time- 
varying  parameters  Is  more  complex  than  the  single-stage 
analysis  described  above,  the  basic  principles  and  conclusions 
remain  valid.  In  particular,  the  generator  performance  Is 
more  or  less  Independent  of  plasma  conductivity  above  200 
mho/cm  but  Is  exponentially  affected  for  conductivities  on 
the  order  of  100  mho/cm  or  less. 

3.2.2  Generator  Design 

Channel  design  la  basically  a problem  of  material 
selection.  It  was  Initially  required  that  the  channel  wall 
material  be  an  electrical  insulator  and  be  sufficiently 
ductile  so  as  not  to  rupture  on  expansion.  Lexan  filled 
these  requirements,  and  was  used  on  Shots  120-1,  -2,  and 
-3.  Later  the  constraint  on  ductility  was  relaxed,  as  the 
entire  channel  was  cast  In  hlgh-denslty  Hydrostone  plaster, 
and  wall  material  was  selected  on  the  basis  of  minimal 
reduction  of  plasma  conductivity  by  ablation  products.  On 
the  basis  of  chemical  equilibrium  calculations  In  the  plasma, 
hydrocarbons  were  shown  to  be  a poor  material,  and  silicon 
dioxide  was  found  to  be  better.  A glass  channel  was  used 
In  120-11. 
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Expansion  of  the  channel  by  plasma  pressures  (^1  kbar) 
was  found  to  be  excessive  In  Shots  120-1  and  120-2,  which 
were  tamped  by  a layer  of  PVC.  Shot  120-3  used  Plaster  of 
Paris  for  tamping,  and  120-^<  used  Hydrostone,  or  high  quality, 
high  density  plaster.  The  plaster  tamping  was  found  to  be 
effective  In  preventing  channel  expansion. 

The  MHD  channel  portion  of  Shots  120-1  and  -2  was 
Intended  solely  to  measure  the  flow  speed  and  gas  conduc- 
tivity (See  Figure  3.7).  An  externally  excited  field  coll 
provided  the  magnetic  field  for  an  open  circuit  voltage 
measurement  (velocity  gage)  and  a loaded  voltage  measure- 
ment (plasma  conductivity).  Since  generator  power  handling 
was  not  a design  consideration,  simple  button  electrodes 
were  used. 

Shot  120-3  was  designed  not  only  to  measure  conductivity 
but  also  to  use  the  generator  to  excite  a field  coll  and  use 
this  field  In  turn  to  drive  a load  (See  Figure  3.8).  To 
maximize  the  gas  conductivity  all  hydrocarbons  were  elimi- 
nated in  the  driver  section,  and  xenon  gas  was  used  Instead 
of  argon.  Electrodes  were  designed  to  protrude  sufficiently 
Into  the  flow  to  penetrate  the  boundary  layer,  and  were 
streamlined  to  minimize  leading  edge  ablation.  The  electrodes 
were  fabricated  from  copper  blocks  to  minimize  circuit  resis- 
tance. In  anticipation  of  a gas  conductivity  of  approximately 
400  mho/cm  the  load  was  designed  to  match  the  expected  Inter- 
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FIGURE  3.7  SCHEMATIC  OF  M H D CHANNEL  CIRCUITRY,  SHOTS  120-1  AND  120-2 
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Diagnostic  Electrodes 


FIGURE  3.8  SCHEMATIC  OF  M H D CHANNEL  CIRCUITRY,  SHOT  120-3 


electrode  plasma  resistance  of  l.^l  mllllohms  in  order  to 
maximize  the  energy  transfer. 


Shot  120-^1  was  designed  to  obtain  full  self-excltatlon 
and  an  efficient  extraction  of  energy  from  the  flowing  gas. 

To  this  end  It  consisted  of  five  switched,  parallel,  self- 
excited  stages  with  a start-up  field  on  the  first  stage  as 

3 

shown  In  Figure  3.9.  Having  multiple  stages  arranged 
sequentially  along  the  channel  and  connected  electrically  In 
parallel  effectively  Increases  the  generator  pulse  duration 
as  seen  by  the  circuit,  and  helps  to  reduce  the  requirements 
on  plasma  conductivity.  Circuit  switching  achieved  by  split 
electrodes  assures  optimum  generator  loading  for  self- 
excltatlon,  and  prevents  overloading  difficulties  at  early 

E 

1' 

times,  when,  for  example,  one  generator  would  otherwise  be  8 

li 

driving  all  five  field  colls.  Although  the  circuit  equations  j 

‘i 

are  complex,  with  time  varying  plasma  resistance,  feedback 
resistance,  and  switching  of  multiple  stages,  conventional  1 

circuit  analysis  techniques  are  applicable.  j 

Electrode  design  of  120-4  was  modified  In  a further  . -j 

t 

attempt  to  reduce  the  observed  plasma  resistance.  The  | 

i 

electrodes  were  finned  to  penetrate  the  channel  wall  boundary  ' 

layer.  They  were  also  plasma-sprayed  with  tungsten  carbide 

on  the  basis  of  calculations  showing  that  ablation  vapor  of 

tungsten  Is  substantially  more  conductive  (a  factor  of  three) 

than  the  ablation  vapor  of  brass. 
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Split  Electrodes 


FIGURE  3.9  SCHEMATIC  OF  MHD  CHANNEL  CIRCUITRY,  SHOT  120-4 


To  minimize  losses  in  external  circuitry  (and  to  prevent 
excessive  heating)  resistances  were  kept  as  small  as  possible. 
At  the  same  time  rlsetlme  considerations  meant  Inductances 
were  to  be  kept  small  and  skin  effect  problems  had  to  be  dealt 
with.  These  considerations  led  to  the  use  of  wide,  flat 
copper  braid  and  wide,  thin  copper  foil  in  field  colls,  loads, 
and  various  interconnections.  The  start-up  coll  and  other 
field  colls  were  kept  as  small  as  possible,  consistent  with 
having  a sufficiently  uniform  field  In  the  electrode  portion 
of  the  channel,  to  maximize  the  magnetic  field  strength  per 
unit  current  flowing  In  the  coll. 

Examples  of  the  experimental  hardware  for  Shot  120-3  are 
shown  In  Figure  3.10.  Figure  3.10  (a)  shows  the  MHD  channel 
before  it  was  cast  In  plaster,  and  Figure  3.10  (b)  shows  the 
driver  and  channel  of  Shot  120-3  assembled  in  the  field  Just 
prior  to  firing.  Figure  3.11  shows  generator  details  of 
120-4,  including  field  coll.  Interstage  circuitry,  and 
electrode  details.  Figure  3.12  shows  the  MHD  channel  as 
assembled  prior  to  casting  and  as  fired. 

3.2.3  Performance  Data 

Performance  of  the  generator  depends  upon  an  effective 
Interaction  between  plasma  parameters  and  circuit  elements. 

The  resistance  and  Inductance  of  each  circuit  element  was 


measured  In  the  laboratory  prior  to  each  shot,  while  the  gas 


(a)  MHD  Channel  With  Coils  And  Diagnostic  Cables  In  Place  Ready  For 
Potting  In  Plaster 


(b)  Driver  And  Channel  Assembly  On  Shot  Stand  Ready  For  Firing 


FIGURE  3.10  PHOTOGRAPHS  OF  EXPLOSIVE  DRIVER  AND  MHD  CHANNEL,  SHOT  120-3 


velocity-history  In  the  MHD  channel  was  measured  during  each 
shot.  Channel  measurements  also  Included  the  time-dependent 
voltage  and  current  at  each  electrode  station.  Prom  this 
data  and  a knowledge  of  the  equivalent  circuit  of  the  MHD 
generator  circuitry  It  was  always  possible  to  determine  a 
unique  plasma  resistance  as  a function  of  time.  This  self- 
consistency  provided  support  for  our  understanding  of  the 
generator  behavior  and  ability  to  predict  electrical  response 
for  different  generator  and  circuit  parameters.  Measured 
plasma  shock  velocity  and  flow  duration  In  the  driver  were 
within  lOJ  of  predicted  values.  From  these  measurements  and 
the  equation  of  state  the  thermodynamic  state  of  the  plasma 
Is  accurat,*ly  known.  Plasma  velocity  Increase  In  the  con- 
vergent section  was  found  to  be  overestimated  by  preshot 
one-dlmenslonal  calculations,  particularly  when  the  area 
convergence  was  larger  than  a factor  of  five.  This  effect 
was  not  unexpected  and  did  not  cause  problems  In  circuit 
design. 

Pull  excitation,  up  to  about  50  tesla  choke  conditions, 
was  not  achieved  In  this  program.  The  difficulty  Is  clearly 
attributable  to  the  plasma  conductivity,  which  was  found  to 
be  substantially  less  than  theoretical  calculations.  The  low 
conductivity  leads  to  a reduced  rate  of  self-excltatlon  of 
the  magnetic  field,  and  affects  exponentially  the  delivered 
power  or  energy  for  a given  plasma  pulse  duration.  Achieved 
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field  amplification  of  about  10,  from  .37  to  3*6  tesla,  demon- 
strated the  concept  of  a self-excited  generator,  but  failed 
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by  a factor  of  15  to  achieve  choke.  On  a logarithmic  scale, 
which  Is  appropriate  for  the  exponential  amplification  process, 
the  achieved  gain  of  20  dB  was  about  half  the  required  gain 
of  43  dB.  Although  the  power  output  was  Impressive  — about 
200  megawatts  — and  broke  all  world  records  for  MHD  power 
production.  It  was  substantially  less  than  peak  design  levels 
of  about  50  gigawatts.  Power  Increases  with  the  square  of 
the  achieved  magnetic  field,  and  full  power  levels  can  only 
be  reached  at  full  field  excitation. 

Initial  theoretical  conductivity  calculations  were  based 
on  the  theory  of  Spltzer  and  Harm  for  a fully  Ionized  low 
density  plasma  (Reference  3*2),  as  modified  later  by  Rogov 
(Reference  3.3)  and  DeVoto  (Reference  3-4).  Careful  experi- 
mental measurements  of  the  conductivity  of  shocked  argon  In 
a shock  tube  (Reference  3*5)  showed  excellent  agreement  with 
theory  up  to  the  maximum  achievable  shock  conditions  In  that 
facility  (a  temperature  of  14,000®  K and  a conductivity  of 
about  80  mho/cm).  This  accumulated  experimental  evidence 
and  the  agreement  among  theoretical  calculations  Indicated 
that  theoretical  conductivity  should  be  accurate  to  well 
within  30X. 

A literature  survey  conducted  under  different  sponsor- 
ship after  technical  completion  of  this  effort  revealed  a 
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significant  Soviet  research  effort  on  dense  plasmas  at 
conditions  very  similar  to  our  experiments.  The  focal  point 
of  the  Soviet  effort  Is  a major  review  article  by  G.  E. 
Norman  (Reference  3.6).  Subsequent  papers  by  other  authors 
have  communicated  both  experimental  and  theoretical  Investi- 
gations of  dense  plasma  properties.  Cesium  was  used  In  many 
of  these  Investigations  because  It  Is  relatively  easy  to 
achieve  the  desired  plasma  conditions.  Some  work  Is  also 
reported  on  argon  and  xenon  plasmas. 

Average  conductivity  was  determined  from  direct 
electrode-to-electrode  resistance  measurements,  relating 
this  value  to  conductivity  by  means  of  an  electrode  constant 
determined  from  measurements  using  copper  sulfate  soultlon. 
Electrode  resistance  was  obtained  directly  from  voltage  and 
current  measurements  at  the  electrode.  Table  3-2  compares 
the  calculated  and  measured  values  of  conductivity  for  the 
tests  In  this  program. 


TABLE  3.2  COMPARISON  OF  PEAK  CONDUCTIVITY  VALUES 


Shot  Number  120-2 

Theory  30? 

Experiment  39 

13 


120-3 

120-4 

Units 

423 

257 

mho /cm 

75 

82 

mho/cm 

.18 

.32 

Ratio 


, Observed  low  conductivity,  and  the  problems  deriving 

from  this  fact  for  self-excltatlon  circuit  design,  led  to  a 
I progressive  series  of  modifications  In  shot  design.  The 

f most  likely  source  of  problems  Is  flow  contamination  from 

I 

( ablating  wall  materials  In  the  driver,  transition  section  or 

channel.  Extensive  calculations  of  thermochemical  equilibrium 
Indicated  that  wall  material  made  a substantatlve  difference 
In  the  conductivity  of  contaminated  plasma.  In  particular. 

It  was  found  that  the  hydrogen  In  hydrocarbons  poisoned  the 
flow  by  Injecting  a large  number  of  atoms  per  unit  mass 
ablated;  equlpartition  demands  that  energy  be  shared  equally, 
and  hydrogen  atoms  rob  energy  otherwise  available  for  Ioniza- 
tion. Ablated  metal  vapor,  on  the  other  hand,  did  very 
little  to  reduce  conductivity.  The  reduction  In  temperature 
was  more  than  equalized  by  the  Increase  In  Ionization  of  the 
metal  atom.  Conductivity  of  the  driver  gas  could  also  be 
Increased  by  using  xenon  Instead  of  argon. 

The  progressive  changes  In  driver  and  channel  design 
to  minimize  contamination  and  boundary  layer  effects  are 
described  earlier  In  this  section.  It  appears  that  elimina- 
tion of  plastic  (hydrocarbon)  material  In  the  taper  section 
and  changing  from  argon  to  xenon  caused  a factor  of  two 
Improvement  In  observed  conductivity,  whereas  further 
Improvements  relating  to  boundary  effects  at  the  channel 
wall  and  the  electrodes  resulted  In  no  significant  Improve- 
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merit  in  observed  conductivity.  The  likely  cause  of  the 
problem  therefore,  although  the  data  Is  scanty.  Is  plasma 
contamination  probably  occurring  In  the  driver. 

There  were  no  apparent  problems  in  the  electrical  circuit 
or  in  the  performance  of  the  generator  given  the  measured 
values  of  plasma  velocity  and  conductivity.  The  data  can 
be  summarized  by  saying  there  were  no  surprises;  generated 
voltages  and  currents  were  consistent  with  the  measured 
values  of  Inductance  and  resistance  In  the  circuit  components. 

Examples  of  the  measured  voltages  and  currents  In  Shot 
120-3  are  shown  In  Figures  3*13  and  3-l^*  Electrical 
schematic  of  this  shot  is  shown  In  Figure  3.8,  and  assembly 
photographs  are  shown  In  Figure  3«10.  The  load  current 
rises  smoothly  to  77  kamps  in  a manner  fully  consistent  with 
the  circuit  parameters. 

Similar  curves  for  Shot  120-4,  which  was  a five-stage 
self-excltatlon  circuit  (See  Figure  3.9)  were  obscured  by  an 
operator-induced  noise  problem  which  caused  baseline  shifts 
In  the  oscilloscopes  (see  discussion  below).  The  backup 
tape  recorder  channels  provided  useful  absolute  data  on 
Stage  3,  and  the  scope  records  provided  data  relative  to  an 
unknown  baseline  for  the  other  stages.  Reduced  data  on  the 
Stage  3 field  coll  current  Is  shown  In  Figure  3.15.  It  rises 
smoothly  from  zero  to  a peak  of  approximately  .29  megamps, 
corresponding  to  a field  strength  of  3.6  tesla.  Peak  power 
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(a)  Field  Coil  Voltage  Versus  Time 


— »i  h* — 20^sec 


(b)  Field  Coil  Current  Versus  Time 


FIGURE  3.13  MHD  EXCITED  FIELD  COIL  VOLTAGE 
AND  CURRENT  HISTORIES,  SHOT  120-3 
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output  of  stage  3 was  in  excess  of  50  megawatts.  Total 
power  output  from  all  five  stages  Is  conservatively  estimated 
to  be  In  excess  of  200  megawatts. 


Field  amplification  did  take  place  In  Shot  120-^4,  but 
was  Insufficient  to  achieve  full  power  output.  Initial  field 
on  Stage  1 was  approximately  0.37  tesla,  and  peak  achieved 
field  In  later  stages  was  approximately  3*6  tesla.  This 
amplification  factor  was  consistent  with  the  observed  con- 
ductivity, but  was  much  less  than  the  factor  of  135  required 
for  choke.  The  peak  measured  conductivity  was  only  82 
mho/cm,  compared  with  the  theoretical  value  of  257 • As 
shown  in  Figure  2.6,  the  amplification  factor  Is  very  sensi- 
tive to  marginal  conductivity  values.  Conductivity  In  excess 
of  150  mho/cm  would  very  likely  have  achieved  full  choke  field 
strengths  of  about  50  tesla  and  power  levels  of  about  50 
gigawatts. 

3.3  Design  and  Performance  of  the  Diagnostic  Devices 

Each  shot  employed  at  least  one  small  set  of  electrodes 
In  a magnetic  field  to  measure  the  velocity  of  the  gas  as  a 
function  of  time.  The  Faraday  open  circuit  voltage  ( v ) 
developed  across  a conductor  moving  with  velocity  ( u ) 
through  a magnetic  field  ( B ) Is  v * Bub  , where  b Is 
the  distance  between  the  electrodes.  Since  B and  b are 
typically  known  parameters,  the  open-circuit  voltage  provides 


a direct  time-dependent  measurement  of  the  velocity  of  the 


Several  factors  influence  the  actual  performance  of  the 
velocity  gage.  Electrode  motion  must  be  small  during  the 
gas  pulse,  otherwise  a proper  b(t)  would  have  to  be  Included 
in  the  analysis.  Fringing  fields  from  nearby  magnetic  field 
sources  must  be  Included  in  the  total  B , and  in  general 
these  might  Include  time-varying  fields. 

The  leads  must  be  Isolated  from  ground  or  at  least 
Isolated  from  other  electrodes  in  contact  with  the  gas  to 
avoid  problems  of  cross-coupling.  This  was  accomplished  by 
using  a wideband  isolation  transformer  between  the  velocity 
gage  and  the  grounded  oscilloscope.  An  example  of  a velocity 
gage  record  is  shown  in  Figure  3.5. 

Examples  of  main  electrode  voltage  records  are  shown  in 
Figures  3.13(a)  and  3.1^(a).  As  with  the  velocity  gage  the 
leads  of  these  had  to  be  Isolated  from  ground  via  wideband 
isolation  transformers,  since  the  main  electrodes  are  like- 
wise in  contact  with  the  gas.  The  only  other  concern  in 
making  these  measurements  was  to  arrange  the  leads  so  as  to 
minimize  inductive  pickup,  since  they  were  generally  in 
locations  of  large  amplitude,  time-varying  magnetic  fields. 

Examples  of  current  measurements  are  shown  in  Figure 
3.13(b)  and  Figure  3.1^(b).  These  measurements  were  made  by 
using  solenoidal  search  colls  closely  coupled  to  current 
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carrying  elements.  The  search  coll  outputs  were  then  inte- 
grated to  provide  signals  proportional  to  circuit  currents 
as  shown.  The  search  colls  were  calibrated  In  the  laboratory 
at  frequencies  corresponding  to  the  MHD  circuit  current  rise- 
times  ('^100  kHz).  The  essence  of  the  calibration  was  to 
measure  the  mutual  Inductance  between  the  search  coll  and 
the  current  carrying  element.  At  the  frequencies  of  Interest 
there  are  negligible  capacitive  effects,  but  care  had  to  be 
taken  to  minimize  stray  inductive  pickup. 

A problem  occurred  in  Shot  120-4  which  is  worthy  of  note 
because  It  affected  the  data.  Shot  120-4  had  numerous 
diagnostic  channels  because  of  the  complexity  of  the  elec- 
trical circuit;  there  were  five  stages,  of  which  three  (the 
first,  third,  and  fifth)  weremonitored  for  field  coll  current, 
electrode  current,  and  electrode  voltage.  In  addition  there 
were  two  velocity  gages.  All  these  components  were  cali- 
brated and  checked  out  In  the  laboratory.  Coupling  problems 
occurred  at  the  test  site,  which  were  finally  traced  down  to 
some  combination  of  ground  shift  in  the  capacitor  power 
supply  for  the  first  stage  saddle  coil,  and  high  impedance 
stray  coupling  through  the  Hydrostone  casting  plaster.  It 
was  found  that  the  problem  could  be  eliminated  by  floating 
the  power  supply  ground.  The  technician  who  fired  the  shot 
Inadvertently  reconnected  the  ground  by  simultaneously 
touching  the  safety  Interlock  key  and  the  power  supply  case. 
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Resulting  base  shifts  affected  the  oscilloscopes  and  most  of 
the  tapedeck  data.  The  problem  could  be  eliminated  In  future 
shots  by  a firing  procedure  change,  a modification  of  the 
power  supply,  or  a careful  redesign  of  the  circuit  Insulation 
In  the  Hydrostone  casting. 
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4.  CONCLUSIONS  AND  RECOMMENDATIONS 
4.1  Accomplishments 


I 


• Designed  and  successfully  tested  explosive  driver 
system.  Conversion  efficiency  from  chemical 
energy  of  explosive  to  delivered  plasma  energy 
was  10$.  Largest  shot  fired  was  120-4,  which 
delivered  12  m.egajoules  plasma  energy  in  approx- 
imately 100  microseconds,  for  a total  thermal 
power  of  120  gigawatts. 

i 

• Designed  and  successfully  tested  fast  self- 
excited  MHD  field  coll  and  circuit  assembly, 
including  multiple  stage  and  switched  electrodes. 

Field  coil  Inductance  for  2-lnch  channel  was 
measured  as  70  nanohenries. 

• Successfully  produced  MHD  power  at  levels  con- 
sistent with  magnetic  field.  Peak  power  was  in 
excess  of  200  megawatts  with  a 3.6  tesla  field 
in  Shot  120-4. 


• Successfully  demonstrated  self-excitation 

consistent  with  the  observed  low  values  of 

i 

electrical  conductivity.  Observed  field 

r 

E 

I amplification  on  Shot  120-4  was  a factor  of 

i 
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approximately  10,  from  .37  to  3.6  tesla. 


• Successfully  developed  and  tested  MHD  measure- 
ment techniques  for  monitoring  high  performance 
energy  flow  under  explosive  field  test  conditions. 
Plasma  velocity  gage  and  electrode  current 
monitors  worked  exceptionally  well  in  all  shots 
but  120-^,  which  had  an  inadvertent  ground  loop 
problem  arising  from  operation  procedures  at  the 
test  site. 

A. 2 Problems 

• Observed  average  conductivity  of  the  plasma  as 
determined  from  electrode  measurements  was  sub- 
stantially smaller  than  theory  (approximately  a 
factor  of  three).  This  low  value,  occurring  In  an 
exponential  term,  prevented  achieving  full  field 
amplification  and  full  power  generation  capa- 
bility. Approximately  1J3  dB  of  gain  (50  tesla) 
was  required  for  full  power,  and  20  dB  (3*6 
tesla)  was  actually  achieved.  An  Increase  of 

50>t  In  experimental  conductivity  would  have 
assured  full  gain  conditions. 


• Experimental  design  modifications  Improved  the 
measured  conductivity,  but  did  not  cure  the 
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problem.  Major  Improvements  followed  changing 
driver  gas  from  argon  to  xenon,  and  substituting 
an  aluminum  convergent  nozzle  for  a PVC  plastic 
nozzle.  Modifications  directed  toward  minimizing 
boundary  layer  effects  — changing  channel  wall 
material  and  electrode  material  and  geometry  — 
did  not  appear  to  Improve  conductivity  signifi- 
cantly. A possible  source  of  the  problem  Is  con- 
tamination of  the  plasma  In  the  explosive  driver 
by  ablated  or  Jetted  wall  material. 

• In  general,  engineering  design  and  fabrication 
of  circuit  components  proved  to  be  considerably 
more  complex  than  Initially  contemplated. 
Calibration  of  the  field  coll  assemblies  and 
MHD  Instrumentation  was  also  more  complex  and 
time-consuming  than  expected. 

Conclusions 

• The  principle  of  self-excited  high-energy,  hlgh- 
power  MHD  electric  power  generation  In  explosive 
tests  scalable  to  large  sizes  has  been  demon- 
strated In  this  program,  but  expected  and  required 
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levels  of  self-excitation  field  gain  for  full 
output  were  not  achieved. 

• In  order  to  achieve  full  power  output  from  a 
self-excited  MHD  generator  the  conductivity 
problem  must  be  solved. 

• Externally-excited  explosive  MHD  is  a reasonable, 
practical  means  for  achieving  extremely  high  plasma 
energy  densities  for  DNA  simulation  requirements 

in  present  and  future  programs. 

^ ^ Recommendat ions 

• The  conductivity  problem  can  be  resolved  in 
either  of  two  different  ways:  identify  and 
rectify  the  plasma  contamination  causing  lower 
than  expected  conductivity,  or  go  to  a con- 
vergent driver  capable  of  delivering  plasma 

at  higher  velocities.  Increased  MHD  feedback 
Impedance  at  high  velocities  reduces  the 
requirements  on  conductivity.  For  DNA  appli- 
cations it  is  recommended  that  the  latter 
convergent  driver  approach  be  pursued.  This 
choice  has  the  advantage  of  operating  at 
Increased  voltage  levels  and  correspondingly 
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more  rapid  self-excltatlon.  It  Is  recommended 
that  a plasma  source  be  developed  with  plasma 
properties  specifically  selected  for  a self- 
exited  MHD  device. 

• Having  demonstrated  a satisfactory  plasma  source 
In  terms  of  flow  velocity,  plasma  conductivity 
and  duration,  a further  effort  Is  recommended 
for  tailoring  both  the  gasdynamlcs  and  elec- 
trical generator  output  to  satisfy  particular 
DNA  simulation  requirements.  This  effort 
consists  of  adapting  the  basic  MHD  technology 
to  meet  the  needs  of  a particular  application. 
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